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Abstract 

We investigate the geometrical features of one-dimensional wave propagation, 

whose dynamics is described by the (2+l)-dimensional Lorentz group. We find 

many interesting geometrical ingredients such as spinorlike behavior of wave 

amplitudes, gauge transformations, Bloch-type equations, and Lorentz-group 

Berry phases. We also propose an optical experiment to verify these effects. 
PACS numbers: 03.65.Bz, 03.40.Kf, 02.40.-k 
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I. INTRODUCTION 



The Berry phases produced by rotations have been most extensively studied for various 
physical systems. Another interesting class consists of Berry phases produced by Lorentz 
transformations. Chiao and Jordan showed that Lorentz-group Berry phases could be 
observed as a change of phase of the electromagnetic field in squeezed states. Successive 
squeezing operations for either light or microwaves by degenerate parametric amplifiers 
induce a phase shift that can be interpreted geometrically. These kinds of experiments have 
drawn considerable interest because they can be viewed as the quantum optical realization 
of the Lorentz kinematics. 

Another optical manifestation of the Lorentz-group Berry's phase was proposed and 
experimentally verified [^. It is shown that when one cycles light through a sequence of 
partially polarized states, the light acquires a geometrical phase. 

It is well known that the wave propagation in one dimension (ID) can be described by 
the Lorentz-group dynamics, therefore, one can expect Lorentz-group Berry phases in wave 
systems. In this paper, we present Berry phases and related geometrical features seen in the 
ID wave propagation in inhomogeneous media. 

In Sec. 0, from the ID wave equation with quasi-periodic potential we derive a 
Schrodinger-like equation, which explicitly exhibits the Lorentz dynamics via the gener- 
ators of SL(2,R). The dynamical variables of the equation are the two quadrature wave 
amplitudes, which form a two-column state vector. 

In Sec. a SL(2,R)-Bloch equation is derived from the Schrodinger-like equation. The 
equation is very similar to the conventional SU(2)-Bloch equation for the spin-1/2 except 
that the SL(2,R)-Bloch vector moves on a conical surface not on a sphere. The similar- 
ity could be very useful for the geometrical understanding of Lorentz dynamics, especially 
because the SU(2)-Bloch model has been successfully employed in various fields such as nu- 
clear magnetic resonances, nonlinear optics, polarization optics, interferometry, and atomic 
physics, etc. 
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In Sec. 1^ we discuss two-valuedness of the state vectors. A Bloch vector represents a 
pair of state vectors with opposite signs. When a Bloch vector circles around the apex of a 
cone, the corresponding state vector changes its sign. It is shown that the two-valuedness 
comes from a symmetry of the wave equation with respect to the translation of half wave- 
length. 

Extending the above geometrical situation, we introduce in Sec. ^ a local gauge trans- 
formation associated with the deformation of space coordinate. 

So far, we have assumed that the state vector is real because the potential function in 
the wave equation is real. It turns out, however, the extension to complex wave functions is 
very convenient, because a complex solution could represent two linearly independent real 
solutions at once. In Sec. 0, we derive a Bloch equation for the case of complex amplitudes. 
In this case, a Bloch vector corresponds to a set of state vectors that differ only in phase 
factors. The Bloch vector moves on the surface of hyperboloid. This is the arena where the 
Berry phase plays a role. In Sec. |VI1| , we define the Berry phase for a cyclic evolution of the 
system. 

Finally, we propose a feasible optical experiment where we can observe the Berry phase 
for ID wave propagation. 

II. ID WAVE EQUATION AND SL(2,R) 

Let us start with a one-dimensional wave equation (the Helmholtz equation) for a 
monochromatic wave \l/(a;)e~^'^* 

r a2 ^ 

+ k^[l + t{x)])m{x)=Q, (1) 



where k = k{uj) is the wave number in free (e = 0) space. 

We assume the spatial variation of the refractive index or the potential e(x) has the 
following form: 

e{x) = €q{x) + 2ec(x) cos 2kx + 2es(x) sin 2kx. (2) 



The spatial modulation at wave number 2k strongly modifies the propagation of the wave 
with wave number k owing to the Bragg effect. The DC component eo also has substantial 
influences through the Bragg resonance condition. 

We represent the wave \E'(a;) with mean wave number k as 

"^{x) = A{x)coskx + B{x)smkx, (3) 

where A{x) and B{x) are slowly varying envelope functions. These envelopes are normally 
considered to be real functions, since e{x) in Eq. (|1]) is real. However, the extension to 
complex functions is sometimes very useful as will be discussed in Sec. [V^. 
Substitution of Eqs. (0) and (^) into Eq. (|I]) gives 
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which can be rewritten in a Schrodinger-like form 

ik-'^\^P) = -e{x)■K\^P). (5) 

with the column vector lip) = {A,B)^. Here we introduced e(x) = (— es,ec,eo) = (ei,e2,e3), 
and K = {Ki, K2, J3), with 

K^ = 1^3/2, K2 = icTi/2, J3 = (T2/2, (6) 

where ai{i = 1, 2, 3) are the Pauli's spin matrices. 

We have assumed that A{x) and B{x) do not change appreciably over distances of the or- 
der of k~^ ( slowly- varying-envelope approximation) and neglected the terms with wavelength 
±3k (secular-term approximation). 

The matrices @ are the generators of SL(2,R), the two dimensional real unimodular 
group, and satisfy the commutation relations 

[K,, K2] = -iJ3, [K2, J3] = iKi, [J3, K,] = iK,. (7) 

The {A, B) plane is squeezed in one direction and stretched in the orthogonal direction by 
the generators Ki or K2, and rotated by J3. 
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The group SL(2,R) has the close connection (locally isomorphic) to the (2+1)- 
dimensional Lorentz group S0(2,l) as does the group SU(1,1). Instead of Eq. (|^), if we 
use a representation 

^(x) = A{x)e-''''^ + S(a;)e''=^ (8) 

then we have the SU(l,l)-Schrodinger equation with K[ = iai/2, K2 = — ia2/2, J3 = (T3/2, 
which satisfy the same commutation relations as Eq. (|^). In the following discussion, we 
mostly use SL(2,R) rather than SU(1,1). 

III. SL(2,R)-BLOCH EQUATION REAL AMPLITUDE CASE 

Equation (|^) has a form similar to the Schrodinger equation for the spin-1/2 [J) in 
time varying magnetic fields B{t); ih{d/dt)\ip) = —B{t) ■ J\ip). The underlying group 
for the spin-1/2 dynamics is SU(2), which is the covering group of the three dimensional 
rotation group, S0(3). The Bloch equation derived from the spin-1/2 Schrodinger equation 
apparently reflects the structure of those groups; the Bloch vector moves on a sphere, and 
the Berry phase for spin-1/2 can be related to the surface area of the sphere. 

In order to see the geometrical structure of SL(2,R)-dynamics, let us derive the corre- 
sponding Bloch equation. (A more detailed derivation will be given in Sec. [V^. See also Ref. 

§•) 

First we define the 2x2 density matrix p as 

p = 2i|7/>)(^|J3 = is-K, (9) 
where represents the row vector (A, 5), and s = (si, 32,53) the Bloch vector. We note, 

s = {2AB,-A'^ + B^A^ + B^). (10) 
The equation of motion for p can be obtained from Eq. (^) as 

ik~'^=[-e.K,p]. (11) 



Using the commutation relations (|^), we have the equation of motion for s 

fc-i^ = ex s, (12) 
dt ^ ^ 

where the SU(1,1) vector product defined as 

a X b = (a2&3 - CLsh, ash - aibs, -0162 + a2&i)- (13) 

for two vectors, a = (01,02,03) and b = (6i,62,&3)- We also define the SU(1,1) scalar 
product 

b = 0161 + 0262 - ash- (14) 

Then from Eq. (p!OD, we have 

s ^ s = 0, S3 > 0, (15) 
which means the Bloch vector s moves on the surface of the (upper) cone depicted in Fig. 

When e{x) is constant, trajectories for s are very simple; the vector s lies in a plane that 
is orthogonal to the vector e = (ei, €2, —63), because 

ds ~ ds _ ^ 

dx dx 

A trajectory, which is the intersection of a plane and the cone could be an ellipse, a 
parabola, or a hyperbola according as e| is greater than, equal to, or less than e{ + e^- 
The elliptical, bound trajectories correspond to the propagating waves in the conduction 
bands and the hyperbolic, unbound trajectories correspond to the evanescent waves in the 
forbidden bands. 

Owing to the geometrical simplicity and the similarity to the conventional Bloch equa- 
tion, the SL(2,R)-Bloch equation could be very useful. 
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IV. TWO-VALUEDNESS 



In this section, we study a subtler geometrical feature of SL(2,R)-Blocli equation. We 
note the correspondence between {A, B) and s given by Eq. (|10|) is not one to one but two to 
one; a single vector s represents two state vectors with opposite signs: {A, B) and {—A^ ~B). 

At first sight, these two state vectors may seem to represent two distinct physical sit- 
uations. Within the framework of the present approximation (slowly varying envelopes), 
however, the two situations are virtually identical. Figure ^ represents the local pictures of 
waves with opposite polarities (real lines). A potential with period ~ ttA;^^ is also shown 
(dashed lines). We note the dispositions of the waves relative to the potential are the same 
and, therefore, the waves would evolve identically. 

This degeneracy comes from the fact that the potential @) is invariant under the trans- 
lation X X + nirk^^ {n: integer), as far as the variation of e{x) over the distance mTk~^ 
can be neglected; 

e(a; + nTrA;^^) ~ e(x). (17) 

On the other hand, from Eq. (^, we see that the translation induces {A, B) —>■ (— l)"(y4, B). 

In order to see the geometrical structure of the two-valued representation, we can make 
a cone from the {A, B) plane as shown in Fig. |^. The (A, B) plane is cut along a line from 
the origin {+A axis in this example) and is wrapped to make a twofold cone with apex angle 
of 60 degrees. When the vector s on the cone encircles a closed curve around the apex, we 
see that the corresponding {A,B) gains a factor (—1)"", where n is the winding number of 
the curve. 

Even though, two vectors (A, B) and {—A, —B) represent (almost) the same physical 
situation, we can conveniently detect the minus sign by interferometric methods, as that for 
spin 1/2 rotated by 27r 0. (The two cases have good analogy, but it is only for the latter 
case that interference experiments are absolutely required and even in principle there is no 
other way to detect the minus sign.) 



7 



Figure ^ shows an example of interferometry. We have two trajectories both of which 
start from {A,B) = (1,0) at x = 0. One evolves with a constant e and the other with 
— e. After some evolution, they reach conjugate points in the {A, B) plane. Looking at the 
corresponding wave forms [Fig. §(b)], we notice that the phase shifts accumulate to make 
signs opposite. 

This situation reminds us of the sign change of a spinor rotated by 27r, the Aharonov- 
Bohm effect with an infinite solenoid, the sign change around the degeneracy of the eigen- 
states of real-Hamiltonian systems P,^, and the (relative) configuration space of two iden- 



tical particles in two-dimensional space ||T0|. In terms of the second analogy, the topological 
magnetic flux of vr^c/e is required to account for the sign change or the phase shift of vr. 
The amount of flux is consistent with the (singular) curvature at the cone apex or the apex 
angle. 

Mathematically speaking, SL(2,R) is a two-valued representation of the (2+1)- 
dimensional Lorentz group, SO (2,1), as SU(2) is for the three-dimensional rotation group, 
SO (3). It is very interesting that the two- valued representation which is believed peculiar 
to the quantum regime manifests itself in the classical context. 

It should be noted that the cone introduced here and that for the Bloch vector discussed 
in the previous section have different apex angles and different parametrization: 

S = (Si, S2, V3S3)/2y/s^, (18) 

where s is a vector in Fig. |^. 

V. GAUGE TRANSFORMATION 

In order to generalize the discussion in the previous section, we introduce a (local) gauge 
transformation. We conveniently use the SU(1,1) amplitudes [see Eq. (^)] A and £: 

A={A + iB)/2, S = e, + ie,. (19) 

With this notation, the Schrodinger equation (^ can be represented in a scalar form as 



2ik'^^A = £A* + eoA, (20) 
ax 

where A* stands for the complex conjugate of A. 

Here we consider a space coordinate x' which is shghtly deviated from x as 

x' = x-^(x). (21) 

We assume the deviation ^ is small enough to assure A{x — ~ A{x), £{x — ^) ~ 
and ^(x — ^) ~ ^(x) for any x. In the x' coordinate, A' and £' are defined as 

a: = Ae-"'^, £' = Se-^"'^ (22) 

respectively, and the equation of motion transforms as 

1 

A'. (23) 



2ik-^-^A' = £'A'* + 
dx' 



° ~'~ dx' 



Here we have an extra term 2d^/dx' which accounts for the gauge transformation. This 
term can be interpreted as a gauge field. 

In Eqs. (I) and (|^), in order to define e. A, and B, we needed a long yardstick accurately 
graduated in 27ik~^. However, with the help of the above equations, we can compare the 
wave amplitudes for two observers whose rulers are not necessarily accurate in a long span. 

We note that the discussion in the previous section can be reproduced by considering the 
case of constant deviation C,{x) = mTk~^, which yields A' = —A, S' = S, and d^/dx' = 0. 
We consider two coordinate systems x and x' shifted by nirk"^, and a third coordinate 
system x" = x + ^{x), with ^{x) = for x < Xi and ,^(x) = nirk^^ for x > X2; in the interval 
of Xi < X < X2, x" smoothly connects x and x'. The integration of the gauge field 

/ = r 2^^dx = 2nnk-\ (24) 
J XI dx 

which is independent of the local behavior of ^, amount to the sign factor exp{—inkl/2) = 
(-1)" between {A, B) and {A',B'). 
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VI. SL(2,R) BLOCH EQUATION COMPLEX AMPLITUDE CASE 



From now on let us consider \%p) = {a, (3)'^ as complex amplitudes. Before proceeding, 
physical meaning of complex amplitudes should be clarified. At first sight of Eq. ([l|), which 
has real coefficients, no complex amplitudes seems required. In fact, if we decompose \ip) 
into real and imaginary parts 

\^p) = \^')+i\r), (25) 

then each of and lip") is a real solution to Eq. (|^). Therefore, a complex solution could 
represent two independent real solutions. The independence is assured by a'P" — a" [3' ^ 0, 
where \'ip') = {a', (3')"^ and = {a", (3")"^. A complex solution is convenient in the sense 
that any real solutions can be represented as Re(z|?/')) with a complex number z. 
Starting from the SL(2, R)-Schrodinger equation with complex 

ik--'-^\^) = N{xm, (26) 

N{X) = eiKi + £2^2 + 63^3, (27) 



let us derive the corresponding Bloch-like equation. The adjoint equation |TT[ is 



ik-'^\^) = N^{xM, (28) 
N\x) = -e,K,-e2K2 + esJ3, (29) 



or equivalently. 



ik-'-^{^\ = {^\N{x). (30) 



In our case, N and iV^ are related via 

= RNR~\ (31) 

where R = ~R-^ = expiinJ^) = 2iJ^. With use of this relation, we can rewrite the adjoint 
equation (|28| ) as 
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ik-^-^Rlcp) = N(x)R\ip). 
dx 



(32) 



and find that R\ip) and obey the same equation. Hence, if R\ip{0)) = \ip{0)), then 
R\ip{x)) = li'ix)) for any x. 

Now we introduce the density operator as 



p=\i:){^\ = \ij){ij\R, 



(33) 



which follows the evolution equation 



The matrix representation of p is given as follows: 

-aP* lap 

It can be parametrized as 



(34) 



(35) 



p = i(so//2 + siKi + S2K2 + S3J3), 



(36) 



where 



So = {a*P — a/3*)/i, 
si = a*P + aP*, 
S2 = -|«|' + |/9|', 



S3 = kr+i/^r. 



(37) 



Substitution of p7| ) and (|36[) into (|3^ ) gives the equation of motion for s = (si, S2, S3): 



(38) 



, -ids 

k — = e X s, 
dx 



which resembles to the Bloch equation, except for the use of the SU(1,1) vector product. 
We note that Sq is a constant of motion and that the relation 
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s ~ s = -So (39) 

holds. 

In the case of sq = 0, the vector s moves on the surface of a cone in the (si, S2, S3) space, 
as in Sec. PTT| ; sq = imphes real amplitudes because a and (3 share the same phase, which 
can be eliminated. 

In the case of Sq 7^ 0, s moves on the surface of a two-sheet hyperboloid, as shown in 
Fig. ^. Without loss of generality, we normalize the state vectors as Sq = 1. 



VII. GEOMETRICAL PHASE 

In the case of real amplitudes, a density matrix represents two states. On the other 
hand, in the case of complex amplitudes, a density matrix p represents an infinite number 
of state vectors. For example, the state \ip') = e~^'^^"\ip), \(p') = e^'^^''\(p) derived from \ip) 
and \ip) with Kg = s ■ K yield the same density matrix; 

p'=|V^')(¥'1=e-'^^1^)(¥^|e^'^^^=P, (40) 

since Kg and p = \il'){'p\ = K^of^ + Kg) commute. In other words, a density matrix or a 
Bloch vector corresponds to an equivalence class of states: |e~"^^^|^/') | — 27r < < 27r|, or 
a ray. 

Now we have a projective structure in which we can introduce the Aharonov-Anandan 
connection [|l^, thereby we can derive the dynamical phase and the geometrical phase 



associated with a cycle of evolution. 

The equivalence class introduced above can be represented more conveniently by using 
a decomposition of the group element of SL(2,R): 

An arbitrary state \ip) can be obtained from \ipo) = 2^-'^/^(l, i)'^, which belongs to a ray 
s = (0,0,1) as 
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(42) 



where \(Pq) = R ^\ipo)- The density matrix for this state is 



-iliJs -iuKi 



(43) 



which does not contain the parameter 0. 

Now we see that all the states in a ray can be derived by applying a group action to the 
states = e~^'^'^'-^\ipo) , whose real and imaginary parts are given as 



— sin 0/2 
COS0/2 



(44) 



cos 0/2 
sin 0/2 

Thus the element in a ray can be parametrized by 0. 

Reflecting the two-valuedness discussed in Sec. fV], the states with and + 27r in a ray 
have opposite signs. 

Let us consider a cyclic evolution in which a state |'0(O)) evolves and returns to a state 
l^PiL)) = exp (i$)|'0(O)), which belongs to the same ray, i.e., p(0) = p{L). The phase 
difference $ is composed of a dynamical phase 0d and a geometrical phase 0g. 

The geometrical phase 0g or the Berry phase for this system is given as [^|13| 



Js 



(45) 



where S is the area enclosed by the path. The integrand is an invariant two form under 
the group action. It should be noted that this two form is different from that for the 
surface curvature of the hyperboloid, (2^3 — l)^^dsi A ds2, or that for the surface area of the 
hyperboloid, (2 - sfY^'^dsi A dss. In the SU(2) case, all of the three quantities have the 
same form incidentally. 



VIII. CONCLUSIONS 



In conclusion, we have explored geometrical features of ID wave propagation. Underlying 
group of this problem is the (2+l)-dimensional Lorentz group. The geometrical approach 
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would be very helpful for intuitive understanding of various wave phenomena. It is also 
applicable to other systems which are governed by the Lorentz dynamics. 

Finally we would like to propose an experiment to observe the Berry phase and other 
geometrical phenomena for the ID wave propagation. In Fig. ^, an optical system is pre- 
sented. We have a medium whose distribution of refractive index e(x) can be controlled 
externally. It is possible to write a holographic grating in optical fibers or in some nonlin- 
ear media. The light propagated through the medium is reflected back with perfect mirror 
M . Thus we can prepare a state corresponding to s(0) = (0, 0, 1) on the mirror side 
{x = 0). By adjusting the mirror position, we can "tune" the phase 0o of the initial state as 
li'iO)) = (cos0o/2, sin0o/2)''" and hence we can map out the complex lip) via Eq. (Q). If we 
set up e(x) so that the evolution be cyclic, then the state on the other side (x = L) should 
be l^PiL)) = (cos02,/2, sin 0^/2)'^. For a cyclic evolution, $ = 0^. — 0o is independent of 0o- 
With the Mach-Zehnder interferometer we can measure the difference $ = 0g + 0d- 
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FIGURES 

FIG. 1. Bloch vector s (for real amplitudes) on the conic surface: sf + ~ ■S3 = 

FIG. 2. Waves with opposite polarities (real lines) have the same phase relationship with the 
half-periodic potential (dashed line). 

FIG. 3. In order to identify {A, B) and {—A, —B), which represent physically the same state, 
we make a twofold cone from the (^4, B) plane by making a cut from the origin and wrapping. 

FIG. 4. Interferometric detection of two-valuedness of wave amplitudes, (a) Two trajectories 
starting off from {A,B) = (1,0) at x = evolved by e = 7r(— 2, 1, 3)/400 (real line), and by — e 
(dashed line), respectively, and reach conjugate points at x = lOOA;"^. (b) Corresponding wave 
forms. 

FIG. 5. Bloch vector s (for complex amplitudes) on the hyperboloid: + — s"^ = —1 
FIG. 6. Proposed experimental setup for measuring the Berry phase for a ID wave system. 
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